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Gravitational wave interferometers are complex instruments, requiring years of commissioning 
to achieve the required sensitivities for the detection of gravitational waves, of order < 10 -21 in 
dimensionless detector strain, in the tens of Hz to several kHz frequency band. Investigations carried 
out by the GEO 600 detector characterisation group have shown that detector characterisation 
techniques are useful when planning for commissioning work. At the time of writing, GEO 600 is the 
only large scale laser interferometer currently in operation running with a high duty factor, ~ 70%, 
limited chiefly by the time spent commissioning the detector. The number of observable gravitational 
wave sources scales as the product of the volume of space to which the detector is sensitive and 
the observation time, so the goal of commissioning is to improve the detector sensitivity with the 
least possible detector down-time. We demonstrate a method for increasing the number of sources 
observable by such a detector, by assessing the severity of non-astrophysical noise contaminations 
to efficiently guide commissioning. This method will be particularly useful in the early stages and 
during the initial science runs of the aLIGO and adVirgo detectors, as they are brought up to design 
performance. 


I. INTRODUCTION 

GEO 600 [1] is one member of a global network of in¬ 
terferometric gravitational wave (GW) detectors aiming 
to make the first direct observation of GWs along with 
the LIGO [ !], Virgo [3], and KAGRA [4] detectors. Cur¬ 
rent state of the art laser interferometric GW detectors 
are instruments designed to detect GWs between ~ 10 Hz 
and a few kHz. GEO 600 is a German-British detector 
located near Hannover, Germany. The optical set up of 
GEO 600 is a dual recycled optical folded arm Michelson 
interferometer with 600 m long beam lines. At the time 
of writing GEO 600 is the only operating large scale in¬ 
terferometer with high duty cycle. Since 2010 a series of 
upgrades, referred to as the GEO-HF scheme, has been 
under way to improve the sensitivity to GWs at frequen¬ 
cies above a few hundred Hz [[]. The upgrades consti¬ 
tuting the GEO-HF program include increasing the laser 
power and the implementation of squeezed light [ >]. Dur¬ 
ing this period of upgrades GEO 600 aims to maintain an 
average duty cycle of ~ 70%, while the remaining ~ 30% 
of time is mostly spent on commissioning activities. 

The frequency-dependent sensitivity of a GW detec¬ 
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tor is limited by a stationary noise floor overlaid with 
non-stationary transient noise artefacts. These compo¬ 
nents limit the sensitivity to GWs signals in two ways: 
a high noise floor masks GWs in the data (reducing the 
signal-to-noise ratio (SNR) of signals), while transients 
can be confused with GWs. These transients produce 
a substantial background of candidate events that limit 
the confidence of GW identification; therefore, we require 
higher amplitude GWs for confident identification. The 
goal of detector characterisation is to identify causes of 
non-stationary noise and the components of the noise 
floor, in order to reduce both. 

GW detectors are complicated instruments made up 
of dozens of subsystems and control loops. Characteris¬ 
ing the detector behaviour works to aid both analysis of 
collected data and commissioning of the detector. The 
GEO 600 detector characterisation (GEODC) group has 
focused mainly on aiding commissioning work, through 
clear and frequent communication with the commission¬ 
ing team, which is powerfully aided by having GEODC 
members on-site at the detector. The GEODC group 
investigates noise artefacts in order to understand and 
eliminate the physical causes. Information is fed back 
to the on-site commissioning team who carry out inves¬ 
tigations, fix hardware issues, and thereby improve the 
instrument sensitivity. The efforts of the GEODC group 
have played a major role in investigations of instrumental 
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artefacts and the removal of their sources. 

We aim to increase the number of potential GW signals 
that are observable in our data. This number is propor¬ 
tional to both the volume of space containing canonical 
signals we are sensitive to and the amount of time for 
which we are observing. By commissioning the detec¬ 
tor we aim to increase the sensitive range, and therefore 
increase the volume of space we are observing. How¬ 
ever, commissioning interrupts the data taking which re¬ 
duces the time spent observing, this time is referred to 
as “down-time”. Every investigation costs observing time, 
so it is important to select them in order to maximize the 
probability of GW detection. 

In this paper we report on a selection of transient noise 
classes that the GEODC group identified between Febru¬ 
ary 2011 and March 2013. The layout of this paper is as 
follows. Section II defines the diagnosis tools used for 
the detector characterisation investigations, and Section 
III describes the procedure for performing these investi¬ 
gations. Section IV briefly describes the key elements of 
GEO 600 that are referred to in the investigations de¬ 
scribed in this paper. Section V explains each noise 
phenomenon investigated. For each investigation, veto 
and commissioning solutions are analysed, when avail¬ 
able, and we use our diagnostic tools to quantify this 
performance. Section VI summarises the investigation 
results, and how these methods should be useful for fu¬ 
ture commissioning efforts. 


II. DIAGNOSTIC TOOLS 

In order to make informed decisions about which com¬ 
missioning interruptions can increase the probability of 
GW detection, it is important to characterise and assess 
the severity of noise sources that are observed in the out¬ 
put data stream of the detector. In this paper we use a 
variety of the standard characterisation software tools for 
this purpose which we describe in this section, however, 
we have developed here a method for utilising them to 
make astrophysically motivated commissioning decisions 
during scientific observing runs of a GW detector. 

We compare initial periods, both with and without the 
use of any available vetoes to reduce the impact of the 
noise in an astrophysical search, with predicted and/or 
real post-commissioning periods. This before verses after 
approach is how we assess the severity of a noise source. 
This informs us of the impact of any planned or per¬ 
formed commissioning interruption. The “predicted” sce¬ 
nario, in which the effect of a potential element of com¬ 
missioning is evaluated in advance, enables decisions to 
be made before performing any commissioning. We give 
results also for the actual impacts of commissioning im¬ 
provements on GEO 600 as a comparison and reflection 
of whether the right choice was made where possible. For 
this paper, however, all of the predicted scenarios were 
created after the fact to demonstrate this procedure us¬ 
ing our knowledge at the time of the noise source to make 


our prediction. The use of vetoes in the initial periods 
will demonstrate how well we can do at removing the 
noise without commissioning, and in some cases offers an 
alternative solution which allows us to relax the need of 
immediate commissioning. 

Strain noise spectral density: is the classic charac¬ 
terisation of stationary noise levels in the detector 
as a function of frequency. When the noise is non¬ 
stationary, the noise spectral density represents the 
stationary components of the noise plus a particu¬ 
lar average of the transients over the time interval 
of the measurement. The spectrum gives the base 
noise level to which any transient (signal or noise) 
must be compared in order to rank its significance. 
It is sensitive to detector configuration changes or 
egregious transients which are either very loud or 
occur very often. See figure 3 for examples. 

O megamaps: give a multiple quality-fact or, time- 
frequency representation of the detector data (for 
example figure 6). However, unlike a spectrogram 
Omegamaps can simultaneously show information 
about both short-duration, broadband; and long- 
duration, narrow-band structures. The map is nor¬ 
malised by the strain noise spectral density to high¬ 
light transient deviation from it. At GEO 600 we 
found this visualisation method to be almost as 
good as the human ear at picking out transient 
events and also added a quantitative element which 
the ears lacked. This tool is used to give a hint at 
which directions to drive the commissioning inves¬ 
tigations and development of vetoes. 

Event trigger generators: produce lists of transient 
noise events for the detector strain data channel 
and a large number of auxiliary channels. This is 
done by thresholding time-frequency maps similar 
to the Omegamaps. These events are referred to 
as “triggers”. The triggers for the detector strain 
data channel are used as one of the inputs to GW 
searches which attempt to separate out those which 
come from signals and those which come from noise. 
At GEO 600 both the hierarchical algorithm for 
curves and ridges (HACR) [ '] and the Omega- 
pipeline [8-1 ( ] are used to produce triggers. Pa¬ 
rameters of the transient are reported, such as the 
central frequency, bandwidth, peak time, duration, 
and SNR. These trigger generators, as well as the 
related Coherent WaveBurst algorithm used in sec¬ 
tion V C, perform similar analyses tasks using dif¬ 
ferent algorithms, and have been used as convenient 
by the GEODC team primarily based on our famil¬ 
iarity with them, and on the availability of support 
from their authors. In principle, a single trigger 
generator could be used, and in practice we relied 
mostly on the Omega-pipeline for the diagnostic 
tools listed below. 
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Cumulative glitch rate histograms: illustrate the 
cumulative SNR rate distribution of triggers (for 
example figure 5(b)). When used on triggers which 
are known to be “noise” with respect to a particular 
search, these histograms can be used along with 
the length of an observation block to calculate 
the false-alarm probability of a candidate signal 
trigger with a given SNR [11]. They thus serve as 
a characterisation of the non-stationary noise in 
the detector and are a necessary intermediate data 
product of every search algorithm. 


Fixed false-alarm probability range: combines in¬ 
formation from the stationary and non-stationary 
noise characterisations, the strain noise spectral 
density and the cumulative glitch rate histogram 
respectively, into a single number which is a mea¬ 
sure of the farthest distance at which a gravita¬ 
tional wave burst source of a given energy could be 
detected in a detector. See [11] for a full descrip¬ 
tion. 

In general the formulation of the range at a given 
frequency / is [ ] 
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where Eqw is the energy released in GWs; p is an 
SNR detection threshold; and S(f) is the strain 
noise spectral density. 

To find the range for a particular frequency band, 
for example the band in which we expect GWs for 
a particular source, we define an integrated range. 
For a frequency band /i to we define the inte¬ 
grated range [ f] as 
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The fixed false-alarm probability range is cal¬ 
culated using the definition in equation 2 with 
p defined by the SNR corresponding to a pre¬ 
determined false-alarm probability level, <a, in the 
cumulative glitch rate histogram. Here all back¬ 
ground triggers that overlap the frequency range 
fi to /2 are used to calculate the histogram. Using 
an SNR defined in this manner assures that we al¬ 
ways define our detection threshold to correspond 
with the given probability, cq of the detection being 
a false positive. 

In practice we use parameters corresponding to 
two putative searches for gravitational wave bursts: 
those emitted from a nearby gamma-ray burst 
(GRB) or those emitted from a galactic supernova 
(SN). For the GRB case we assume an emission 
of Eqw — 1O _2 M 0 c 2 in the frequency band of 


100 Hz-500 Hz, for the SN case we assume an emis¬ 
sion of Eqw — 10 ~ 8 MqC 2 in the frequency band of 
500 Hz-4 kHz. In both cases we use an SNR thresh¬ 
old corresponding to a false alarm probability of 
^ 10 -3 in a few seconds long search window around 
the neutrino or gamma-ray observation of the SN 
or GRB, with a corresponding false alarm rate of 
5.6 x lO -4 Hz. These two ranges give us the basic 
astrophysical figure-of-merits for the low and high 
frequency sensitivity of GEO 600. 

Spacetime observation rate: is an extension of the 
fixed false-alarm probability range which measures 
the rate at which we are accumulating observed 
spacetime volume for a particular source. This 
is derived from the quantity of most interest to 
our astrophysical searches: the expected number 
of observed sources, A^qw, during a period of ob¬ 
servation of length t. For a uniform distribution of 
sources, H, as in [ ] we have 

JV GW = y (RfBt. (3) 

Expanding the observation time, £, as the product 
of the total time, T, and an operation duty cycle, 
D, we can then calculate the spacetime observation 
rate as 


= (4) 

This is the basic figure-of-merit we will use to com¬ 
pare two configurations of data production in a GW 
detector. It combines information about the sensi¬ 
tivity of the detector through the fixed false-alarm 
probability range, (i2), with the duty cycle of data 
production, D, hindered either by regular detector 
down-times or time-domain vetoes. 

Vetoes: consist in removing from the analysis triggers 
that are in a contaminated time-frequency regions. 
These contaminated time-frequency regions are de¬ 
rived from investigations and auxiliary sensors in 
the detector. We use two types of vetoes: time- 
domain vetoes and frequency-domain vetoes. Ve¬ 
toes remove glitches and hence reduce the p thresh¬ 
old in equation 1, however this comes at a cost that 
depend on the veto type. 

Time-domain vetoes remove triggers with central 
times within specific slices of time. Hence they 
reduce the effective duty cycle D. Frequency- 
domain vetoes remove triggers with central fre¬ 
quency within specific frequency bands. Hence they 
render the detector blind at some frequencies, the 
range R(f) in equation 2 is set to zero in the vetoed 
frequency bands. 

The spacetime observation rate includes both the 
benefits and the cost of using vetoes. 
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These benefit-delay time figures-of-merit are the 
quantities which we will use to drive commissioning 
decisions. If the benefit-delay time is shorter than 
the remaining time for an observing period after a 
commissioning break, then the commissioning will 
increase observable spacetime during the observing 
period. If this is not the case, the loss in observ¬ 
able spacetime due to commissioning will not be 
re-accumulated during the observation period. For 
this reason, it is important to estimate t comm con¬ 
servatively, taking the most pessimistic reasonable 
estimate into account. 


FIG. 1. Cartoon showing the spacetime observation rate (5) 
against time for a noise investigation. As defined in equation 
4, areas under this curve are proportional to the expected 
number of observed sources. The blue area denotes a pe¬ 
riod of detector performance while the noise was present be¬ 
fore the commissioning of an improvement occurs. In green 
we show a period after the noise was removed. The grey 
area under the dashed spacetime observation rate denoted by 
“I” represents the potential observed sources which were lost 
due to the break in operation caused by commissioning. The 
benefit-delay time Tbenefit-deiay, is then defined as the time 
required, while running the detector at the improved space- 
time observation rate after commissioning, to make up for the 
commissioning loss. As the area denoted by “11” represents 
the improvement in spacetime observation rate, Tbenefit-deiay 
is the time which makes Area(II) = Area(I). 

Benefit-delay time and benefit-delay ratio: are 

figures-of-merit which help us determine how long 
we need to operate a detector at an improved 
sensitivity to make up the deficit number of 
observable sources during the period when the 
detector is off-line for commissioning. 

The benefit-delay time, Tb en efit-deiay is depicted in 
figure 1. Here we see that 

™ - S b* comm /r\ 

-^benefit-delay — "V V - (pj 

5 a 5b 

where £ com m is the estimated down-time due to 
commissioning and 5b and 5 a are the spacetime 
observation rates before and after the commission¬ 
ing period respectively. It should be noted that this 
definition assumes that the spacetime observation 
rate is constant enough before the commissioning 
interruption that the deficit number of observable 
sources caused by commissioning is proportional to 

comm • 

It is useful to express the benefit-delay time figure- 
of-merit as the ratio 

^ _ ^benefit-delay 

^benefit-delay — 7 •> vP/ 

komm 

so that the benefit-delay time can be easily calcu¬ 
lated from any commissioning down-time. 


III. PROCEDURE FOR INVESTIGATIONS 

The broad layout of each investigation in section V 
is the same, although details vary in each due to the 
specifics of the situations. This section describes the 
common approach, and how to use the tools described 
in section II to plan detector commissioning. The in¬ 
vestigations presented in this article were analysed retro¬ 
spectively, except for the final case of the 3.5 Hz dither¬ 
squeezing glitches, in section V C, in which a prototypical 
analysis was carried out at the time to make just such a 
commissioning decision. We nonetheless outline a proce¬ 
dure for solving noise problems, as they arise, in a manner 
which is most efficient with respect to the astrophysical 
productivity of a GW detector. 

Two standard techniques are available to mitigate noise 
contamination of the data. The first approach is to re¬ 
move the source of the noise through commissioning the 
detector. The second approach is to excise the data 
affected in the contaminated time-frequency band from 
analysis, referred to as vetoing. Vetoes can be applied af¬ 
ter the fact, and do not incur the substantial additional 
periods of down-time from commissioning. However, ef¬ 
fective vetoes are not always available if the phenomenon 
is severe enough to be endemic or subtle enough to evade 
automatic veto generation. Commissioning removes the 
noise source, meaning that the only lost data is from a 
one-time commissioning effort. Unfortunately, the down¬ 
time required for investigations and commissioning can 
be extensive and difficult to predict, and successful noise 
reduction is not assured. 

Initial identification of a problematic noise source is 
generally through the effect on the spacetime observa¬ 
tion rate. This either results from the effects on the 
strain noise spectral density, the cumulative glitch rate 
histogram, or both, depending on the stationarity of the 
noise source. These scenarios are distinguished by plot¬ 
ting the strain noise spectral density and the cumulative 
glitch rate histogram as in figure 5. Non-stationarity is 
further exposed by Omegamaps as in figure 6, which show 
the transient nature of the noise, providing more infor¬ 
mation to assist the determination of the origin of the 
noise, and also contributing to the possible development 
of any veto. 
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After identification of the cause of the noise, investi¬ 
gations on the data can proceed, initially without inter¬ 
fering with detector operation. During this time plan¬ 
ning for experimental tests and possible commissioning 
solutions are made. These investigations build an under¬ 
standing of the noise extent so as to predict the improve¬ 
ment in spacetime observation rate without the noise, 
and when possible also to produce a veto. Additionally, 
short investigations with the instrument are often nec¬ 
essary to understand the noise well enough to generate 
these elements. 

We can then apply the benefit-delay time formalism, 
defined in section II by equation 5 and equation 6. Here 
we use the predicted spacetime observation rate for the 
“after commissioning” scenario to quantify the expected 
benefit from commissioning the noise source. If a veto 
can be constructed, we then calculate the spacetime ob¬ 
servation rate with the veto applied and substitute the 
“before commissioning” scenario in the benefit-delay time 
analysis with this vetoed scenario. The comparison of the 
spacetime observation rates from these analyses results in 
a measure of how effective any available vetoes can be in 
mitigating the need for commissioning. 

For the rest of this article we work through the gen¬ 
eral method by describing important noise contamina¬ 
tions, and the associated commissioning campaigns in 
GEO 600. For two of the cases in this article, we have 
easily comparable before- and after-commissioning sce¬ 
narios where only one problem was addressed during the 
commissioning down-times. In situations like this, it is 
possible to validate the benefit-delay time analysis, to re¬ 
veal if the commissioning decision was appropriate. For 
the final case, we discuss the use and investigate the cor¬ 
rectness of a prototype version of this process that was 
applied in order to decide against aggressive commission¬ 
ing. 


IV. GEO 600 OVERVIEW 

Before going into the investigations in section V, it 
is important to give a brief description of some of the 
relevant elements of GEO 600. A simple schematic of the 
optical layout of GEO 600 is illustrated in figure 2. Laser 
light at a wavelength of 1064 nm is fed into the vacuum 
system from an amplitude and frequency stabilised laser. 
It first passes through two suspended triangular optical 
cavities that remove higher order spatial modes of the 
laser light and reduce amplitude noise. These cavities 
are called the input mode cleaners. 

The interferometer itself is composed of seven large 
suspended main optics, which form a power and signal re¬ 
cycled Michelson interferometer with 1200 m long folded 
arms. 

At the output of the interferometer, three beam di¬ 
recting optics (the BDO mirrors) are used to direct the 
light into a Faraday isolator and through an output mode 
cleaner (OMC). The purpose of the BDOs is twofold: 


MFn 



FIG. 2. Simple layout of GEO 600, see section IV for details. 
Included are the BDOs, OMC, and squeezer which are rele¬ 
vant for the investigations in section V. 

they form a telescope which reduces the large beam waist 
of the interferometer to match the small waist of the out¬ 
put optics while also serving as alignment actuators. 

The OMC rejects radio frequency side-bands of the 
laser light that are used for interferometer control and 
also rejects higher order spatial modes coming from de¬ 
fects of the main optics. The power in this cleaned light 
exiting the OMC is detected by the photo-diode shown. 
This signal is used to measure changes in the relative arm 
length inside the interferometer, and is referred to as the 
strain data channel. 

The Faraday isolator is required in the procedure used 
for the injection of a squeezed vacuum state into the out¬ 
put port of the interferometer. This squeezed light de¬ 
creases the shot noise at the photo-diode, the dominant 
source of noise in the relative arm length sensing at high 
frequencies (see section VB for details). The sensitivity 
of GEO 600 can be seen in figure 3 where the amplitude 
spectral density of detector strain is plotted for different 
epochs that span the work reported in this paper. The 
amplitude spectral density is just one way of quantifying 
the astrophysical sensitivity of a GW detector. 

V. INVESTIGATIONS 

This section details three investigations into noise 
sources that occurred at GEO 600. This is not an ex¬ 
haustive summary of the work done by either the detec¬ 
tor characterisation or commissioning teams. The data 
used in this paper is not representative of the general 
sensitivity of GEO 600, rather these investigations were 
chosen to display a variety of noise sources requiring dif¬ 
ferent balances of vetoes and commissioning. The “be¬ 
fore” and “after” commissioning times for each investi¬ 
gation were chosen to highlight the problems clearly, so 
ideally only the effects of the noise sources being investi¬ 
gated are present. However, there are many noise sources 
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Effect of Commissioning on GEQ600 Sensitivity 


— January 2011 
June 2011 
February 2012 

— July 2013 



Frequency (Hz) 


FIG. 3. Strain noise spectral density versus frequency for GEO 600 for different epochs referred to in section V. In black is 
depicted the state before the characterisation and commissioning campaign of this paper; in cyan after introducing additional 
seismic isolation stages to the BDOs; in green after additional isolation of the OMC; and in red after removal of other transient 
noises and narrow line features, as well as improvement of squeezer performance. 


whose effects on the strain data channel stream show up 
and disappear unpredictably so this will not always be 
the case. 

A key element for why these investigations were chosen 
is their ability to demonstrate how a close relationship be¬ 
tween a characterisation team and a commissioning team 
can improve the scientific output of a GW detector. Im¬ 
portant insights are obtained by retroactively applying 
the robust analysis outlined in section III to key noise 
investigations, showing quantitatively how detector char¬ 
acterisation can and should contribute significantly to the 
direction of the detector commissioning effort; resulting 
in an overall increase in the number of observable sources 
compared to unguided commissioning. 

The remainder of this sections details our investiga¬ 
tions into different noise sources observed at GEO 600. 
Each investigation discussion describes the noise phe¬ 
nomena, veto and commissioning results, compares their 
performance, and discusses the best solution and other 
lessons learned. 


A. OMC alignment issues 

The three optics that match the main interferometer 
to the OMC are referred to as the BDOs. Originally the 
BDO suspensions were not designed to be compatible 
with an OMC for DC readout, so they were hung as sin¬ 
gle pendulum suspensions from blade springs, designed to 
provide isolation in the vertical direction. These springs 
are mounted on a rigid support within the vacuum cham- 



FIG. 4. The transmission power, P, of the beam exiting the 
OMC as a function of an input beam alignment degree of 
freedom parameter 0 is shown as the black curve. A high 
frequency fluctuation in 0 shown in red (a) is modulated by a 
low frequency variation shown in green (b). This produces a 
modulated high frequency fluctuation in the transmission (c) 
shown in red with the low frequency window shown in green. 


ber and are effectively connected to the ground at the 
frequencies of interest. Due to this single stage design, 
ground motion easily excites resonances of this system, 
coupling into mirror motion and ultimately misalignment 
of the beam impinging on the OMC. 

The transverse motion of the suspension point of the 
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optics hung on pendula is highly suppressed at frequen¬ 
cies above the resonance frequencies of the pendula. At 
frequencies well below the resonance frequencies there is 
no such suppression, and there is amplification near the 
resonances. At GEO 600, for the main interferometer the 
optic’s suspensions have resonance frequencies spanning 
0.5 Hz-2 Hz. At these frequencies there are numerous un¬ 
controlled interferometer degrees of freedom. This leads 
to slow drifts of, for example, the position and angle of 
the output beam of the interferometer. In turn, this mo¬ 
tion misaligns the output beam onto the OMC. 

Any additional high frequency alignment fluctuations 
are modulated by these low frequency variations, produc¬ 
ing a noise source that varies in amplitude as the align¬ 
ment of the OMC changes. This mechanism is shown in 
figure 4. When the OMC is misaligned, a fluctuation en¬ 
tering the OMC will couple into the amplitude transmis¬ 
sion of the OMC, and thus the observed detector strain 
signal, at the frequency of the fluctuation. Conversely, 
when the OMC is well aligned, the entering fluctuation 
will couple more weakly and only to the higher order 
harmonics of the frequency of fluctuation. 

Of particular concern were the violin modes of the steel 
wires suspending the BDOs which spanned frequencies 
from 140 Hz-205 Hz. Resonances of the support struc¬ 
ture from which BDOs are suspended spanned a range of 
frequencies from 80 Hz-600 Hz. Through the mechanism 
described above, these high frequency fluctuations that 
show up in the alignment of the OMC were modulated 
by the large low frequency variation due to the motion 
of the main optics. 

This produced complex noise features in the strain 
data channel mainly seen from 100 Hz-500 Hz; see fig¬ 
ure 5(a). When comparing noise spectra from before 
and after commissioning, we see that at some frequen¬ 
cies there is an increase in the noise by almost a fac¬ 
tor of 100. This noise is also seen very strongly in the 
cumulative glitch rate histogram shown in figure 5(b). 
For SNRs > 20 we see an increase of over a factor 100 
in the glitch rate. More importantly, the SNR thresh¬ 
old at which background triggers happen at a rate of 
5.6 x 10 _4 Hz, which is a factor in the computation of 
the spacetime observation rate, increases by a factor of 
two. 

It should be kept in mind that at the time these effects 
on the strain data channel were observed, the details of 
the source of the noise and coupling mechanism were un¬ 
known. 

Plotting the Omegamap of the strain data channel as 
shown in figure 6 displays the nature of the noise more 
clearly, by highlighting two aspects which add to the in¬ 
formation provided by the characterisations in figure 5. 
The first is that the noise is confined to relatively well- 
defined frequencies which remain constant. The second 
is that the noise is highly variable, at time scales rang¬ 
ing from 0.5-2 s, with occasional times where all resonant 
bands (near 140 Hz, 170 Hz, 200 Hz, and 350 Hz) are si¬ 
multaneously quiet. With the knowledge of the source 


of the noise, it is clear where these characteristics come 
from. The well defined frequencies are jitter frequen¬ 
cies of the BDO suspensions while the transient nature 
comes from the motion of the main optic suspensions of 
GEO 600. Here the quiet times are short periods when 
the OMC passes through the correct alignment and the 
jitter line coupling is suppressed. 

It is already clear from these initial views of the prob¬ 
lem how to move towards designing a veto which could 
be used to mitigate its effects. We use the well-defined 
frequency aspect of the noise to create frequency-domain 
vetoes by notching various frequency bands. The notched 
bands from 74 Hz-110 Hz, 130 Hz-148 Hz, 163 Hz-242 Hz 
and 343 Hz-356 Hz remove 30% of the bandwidth be¬ 
tween 100 Hz-500 Hz. Vetoing these frequency bands re¬ 
moves the noise source from the data, but also makes the 
detector blind to GW signals in these frequency bands. 
When carrying out an analysis or calculating the space- 
time observation rate, we apply these notches to both the 
strain noise spectral density and the triggers. For the 
spectrum, it is assumed that the noise amplitude within 
the notch band is infinite. The triggers are notched by 
removing the glitches with central frequencies within the 
notched bands. The effects of this veto on the spectrum 
and cumulative trigger rate are shown in figure 5(a) and 
5(b) respectively. The SNR threshold is reduced by ap¬ 
proximately a factor of two. Table I displays the space- 
time observation rates for all investigations discussed in 
this article. It shows there that this veto increases the 
spacetime observation rate by a factor of 8. 

For this investigation, there was a straightforward way 
in which a prediction for the effects of commissioning im¬ 
provements could have been made. This noise was seen 
on switching the readout technique from a heterodyne- 
readout to a DC-readout [ ]. Therefore it could be 

assumed that with commissioning it would be possible 
to reduce the noise back down to the level seen with 
GEO 600 configured in the heterodyne-readout scheme. 
Thus we take one of the heterodyne-readout strain noise 
spectra in the frequency band of interest as a predic¬ 
tion of the commissioning improvements. This is shown 
in figure 5(a). We combine this spectrum with the noise 
transients from the before commissioning period with the 
frequency vetoes applied to calculate the predicted range. 
This gives us a predicted burst-like spacetime observation 
rate of 300 x 10 5 kpc 3 as can be seen in table I —a factor 
of 50 improvement over the before-commissioning obser¬ 
vation rate and a factor 7 over the with-vetoes rate. 

With vetoes and a prediction for improvements from 
commissioning, we can carry out a benefit-delay analy¬ 
sis as described in section II. The results are given in 
table II, as well as in figure 7 where we have already ap¬ 
plied the actual amount of time which was needed for 
commissioning the improvements. Here we see that al¬ 
though the veto already provides very substantial im¬ 
provements to the spacetime observation rate, the predic¬ 
tion promises such large reductions in noise that the as- 
trophysical losses from commissioning is made up within 




(b) Cumulative glitch rate 


FIG. 5. BDO suspension upgrade comparison plots, in black we show a period of data before the BDO suspension upgrade while 
a period of data after the upgrade is shown in green. The improvement in the noise spectral density from the commissioning 
can be seen as the reduction of the complex noise features from 100 Hz-450 Hz in (a). The glitch rate in the 100 Hz-500 Hz band 
with SNR > 20 decreased by two orders of magnitude, as shown in (b). Also, the SNR threshold, where louder background 
triggers happen only at a rate lower than 5.6 x 10 -4 Hz (horizontal dashed black line), is reduced by approximately a factor 
of two. In red we display the effects on these noise characterisations of the application of vetoes defined in the text, while 
turquoise shows our prediction of the spectral noise reduction from commissioning. 



FIG. 6. Omegamap showing a period of time before the BDO 
suspension upgrade. Low frequency (0.5 Hz-2 Hz) fluctua¬ 
tions in the alignment of the OMC from motion in the inter¬ 
ferometer couple BDO suspension jitter in the few hundred 
Hz range (i.e. 140 Hz) into the detector strain data channel. 


a fraction of the time taken to do the commissioning. 
Since at this time we knew that GEO 600 would be run¬ 
ning for at least three full years more, commissioning an 
improvement would have to take more than 2.5 years of 
lost operation time for it to be no longer worthwhile. In 
this case it would have been clear that commissioning 
interruptions should proceed. 

From the spectral characteristics of the noise alone, 
however, invasive commissioning investigations were al¬ 
ready deemed necessary before seeing such benefit-delay 
time analyses although it was unclear how much time 
would need to be devoted to find the source of the noise. 
Once the source was identified as described in the be¬ 
ginning of this section, the BDO suspensions were re¬ 
designed to provide roughly 100 times better isolation of 


the ground motion from 50 Hz-1 kHz. This reduced the 
high frequency fluctuations which were being modulated 
to create the large noise structures in figure 5(a). The 
cumulative down-time from the investigations as well as 
commissioning of this upgrade was 20.8 days. The ef¬ 
fect of this commissioning on the spectrum and cumu¬ 
lative trigger rate are shown in figure 5(a) and 5(b) re¬ 
spectively. The transient noise shows improvements very 
similar to the improvements gained by applying the ve¬ 
toes. For the spectrum we also observe large gains com¬ 
pared to the noise level before commissioning. Firstly, 
the noise is reduced over almost the entire band. Sec¬ 
ondly, many of the large noise features are significantly 
reduced: in the best case, near 200 Hz, by almost a fac¬ 
tor of 100. These improvements together contribute to 
a 17-fold increase in the spacetime observation rate over 
the before-commissioning rate and only a 2-fold increase 
for the vetoed-case. 

Using these actual results of the commissioning im¬ 
provements we can again carry out a benefit-delay time 
analysis to see if the commissioning interruptions were 
actually worthwhile. In table II and figure 7 we see that 
although the improvement turned out not to be better 
than the prediction, the deficit in potentially observed 
sources caused by the commissioning is very quickly made 
up by the improvements seen. The benefit-delay ratio 
for the actual improvements, with vetoes applied, is 0.75 
so we would have always recovered from the number-of- 
sources deficit in less time than the total commission¬ 
ing down-time. Through our cost-benefit analysis we see 
that the commissioning interruption was worthwhile. Al¬ 
though this was the assumed outcome simply by looking 
at the spectral noise amplitude, the quantitative nature 
of this analysis justifies this assumption. In the follow¬ 
ing examples, we will see cases where these benefit-delay 
time analyses would have had an impact on the way we 
thought about the commissioning—either because of the 
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Time [days] 

FIG. 7. Graphical display of the spacetime observation vol¬ 
ume for a GW search, S = f S dt around the time of the 
BDO suspension upgrade. The period displayed here starts 
from sometime before the commissioning period until after 
all benefit-delay times have been reached. All spacetime ob¬ 
servation volumes are shifted so they are equal at the be¬ 
ginning of the commissioning period which is designated by 
the grey filled in area. Here all real spacetime observation 
volumes are plotted as solid lines while extrapolated or pre¬ 
dicted spacetime observation volumes are plotted as dashed 
lines. In black we plot quantities which correspond to the 
before-commissioning period without vetoes while red denotes 
those quantities with vetoes. The green lines plot spacetime 
observation volumes during and after commissioning (hypo¬ 
thetical and real). The benefit-delay times are displayed by 
the vertical lines where the hypothetical continued before- 
and veto-spacetime observation volume meets the after- or 
predicted-spacetime observation volume. 


subtlety of the problem or the power of the vetoes. 

We take a look now at a noise source involving the 
same mechanism as the noise caused by the BDO sus¬ 
pensions but with far subtler effects on the strain data. 
This source has resonances of elements within the OMC 
and its suspension. The resultant noise had been previ¬ 
ously buried by the noise caused by the BDO suspensions 
so was only noticed after the BDO suspension upgrade. 
Here we observed resonant structures near 255 Hz and 
280 Hz which can be seen by comparing the before- and 
after-commissioning spectra in figure 8(a). These struc¬ 
tures are not more than a factor 4 above the underlying 
noise floor. In the cumulative glitch rate histogram (fig¬ 
ure 8(b)), triggers with SNR > 10 occur a little less than 
8 times more often when the noise is present than when 
it is not present and the SNR threshold is about 30% 
higher. 

The nature of the noise can be seen in the Omegamap 
of the detector output channel in figure 9. Comparing 
the Omegamaps figure 6 and figure 9, without previous 
knowledge of the mechanisms, we see that these noise 
sources have similar frequency characteristics (both sta¬ 
ble) and similar transient behaviour. From these plots 
one may already guess that the mechanism for coupling 


of motion into the strain data channel for both noise sce¬ 
narios is the same. 

As before, a frequency dependent veto can notch out 
the glitch bands while maintaining sensitivity at other 
frequencies. The veto used here is a notch at 243 Hz- 
293 Hz applied in the same manner as done above. This 
time the veto only removes 13% of the bandwidth be¬ 
tween 100 Hz-500 Hz. The effects of this veto on the 
spectrum and cumulative trigger rate are shown in fig¬ 
ure 8(a) 8(b) respectively. The SNR threshold is only 
reduced by 16%. In table I we see that the spacetime 
observation rate increases by 60% through the use of the 
veto. 

With the hint from the Omegamap that this noise orig¬ 
inates from a mechanism similar to that of the BDO sus¬ 
pension noise, it is not difficult to come up with a predic¬ 
tion on what we can achieve with commissioning. Here, 
since only a factor of 5 is needed, we assume that we 
can suppress these resonant structures to a level where 
they would lie below the underlying noise floor. We then 
take as a predicted spectrum, one where we have sim¬ 
ply interpolated across the notch as shown in figure 8(a). 
Combining this spectrum with the noise transients from 
the before-commissioning period, with the frequency veto 
applied, then allows us to calculate the predicted range. 
This gives us a predicted burst-like spacetime observa¬ 
tion rate of 850 x 10 5 kpc 3 as can be seen in table I— an 
80% improvement over the before-commissioning space- 
time observation rate but only a 13% improvement over 
the with-vetoes rate. 

Carrying out a benefit-delay time analysis using 
the spacetime observation rates from the vetoed- and 
predicted-scenario then yields the results in table II, fig¬ 
ure 10 where we get a benefit-delay ratio of 8. This ra¬ 
tio is much larger than for the noise which called for 
the BDO suspension upgrade, however, still low enough 
that commissioning interruptions for investigations and 
improvements would be beneficial. Even if the commis¬ 
sioning down-time would have taken as long as the BDO 
suspension upgrade down-time, the deficit in the poten¬ 
tially observed GW sources would be made up within a 
half year of observation. However, due to observations of 
a change in the resonant noise structures after an instru¬ 
mental change commissioned for other reasons, a little 
after this problem was observed it was already clear that 
a solution would require much less time. The commis¬ 
sioning interruption was thought to take no more than 
3 days which would bring the benefit-delay time down 
to 3.5 weeks. Given these considerations, it is clear that 
commissioning a solution to this noise should proceed as 
soon as possible. Due to the subtle nature of the problem, 
without the benefit-delay time analysis this conclusion is 
not obvious. 

The serendipitous observation of a change in the two 
resonant structures happened directly after work on the 
OMC optical bench. Thus we decided to improve the iso¬ 
lation of this bench from ground motion as an attempt 
to remove these noise structures. The cumulative down- 
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(b) Cumulative glitch rate 


FIG. 8. OMC isolation upgrade comparison plots, in black we show a period of data before the OMC isolation upgrade while 
a period of data after the upgrade is shown in green. The noise features are well confined to two narrow frequency bands 
around 255 Hz and 280 Hz, where we can easily select a notch veto (red) to remove these features from the spectrum, as can be 
seen in (a). The glitch rates after applying vetoes and after commissioning are similar: both are about a factor 5 lower than 
the before-commissioning level for SNRs > 20 as can be seen in (b). However, the SNR threshold, where louder background 
triggers happen only at a rate lower than 5.6 X 10 -4 Hz (horizontal dashed black line), is reduced only by 20%. The turquoise 
curve in (a) shows our prediction of the spectral noise reduction from commissioning. 



FIG. 9. Omegamap showing a period of time with transient 
noise in the two bands around 255 Hz and 280 Hz. Fluctua¬ 
tions in the alignment of the OMC allow noise entering the 
detector through the seismic isolation of the OMC to be cou¬ 
pled into the strain data channel. 


time from the required investigations and commissioning 
work was 1.9 days. The effect of this commissioning on 
the spectrum and cumulative trigger rate are shown in 
figure 8(a) and 8(b) respectively. The amplitude of the 
resonance structures at 255 Hz and 280 Hz are reduced 
by a factor of almost 3 and 2 respectively, see figure 8(a), 
such that they are beneath the noise floor. We also see 
small improvements in other features in the noise floor 
across the 100 Hz-500 Hz band. Similar but a little more 
effective than the application of the veto, SNR threshold 
reduces by a factor of 22%, as seen in figure 8(b). These 
small improvements lead to a 3-fold increase in the space- 
time observation rate over the before-commissioning rate 
and a 70% increase over the vetoed-case. Comparing 
these improvements to the predicted values in table I 


shows that we have, this time, made conservative predic¬ 
tions. 

Putting the actual spacetime observation rate into the 
benefit-delay time analysis then highlights the fact that 
the commissioning brought better improvements than 
predicted by lowering the benefit-delay ratio from 8 to 
1.5, as shown in table II and graphically in figure 10. 
This shows that the decision to proceed with commis¬ 
sioning was in fact well justified. 

In contrast to the noise-case solved by upgrading the 
BDO suspensions, we show that the decision to commis¬ 
sion an upgrade here would have attracted more weight 
had a benefit-delay time analysis been carried out. In 
the next section we will look at an example for which a 
benefit-delay time analysis gives key insight into the best 
solution for a severe transient noise source. 


B. Squeezer glitches 

GEO 600 is the first large scale interferometric detec¬ 
tor to implement the injection of squeezed light to im¬ 
prove the shot noise limited sensitivity [14], and also the 
first such detector to operate this technique on year time 
scales [ 5] and during science operations [16]. 

To gain an understanding of squeezing, we can pic¬ 
ture the quantum uncertainty by drawing an ellipse in a 
plane spanned by two conjugate variables, which for laser 
light are the phase and amplitude of the electric field as 
shown in figure 11. The minimum area of this ellipse 
is set by the Heisenberg uncertainty principle. Squeezed 
states of light are quantum states that have a reduced 
uncertainty in one of the field quadratures compared to 
a coherent state, and an increased uncertainty in the con¬ 
jugate quadrature. A squeezed state of light is injected 
into the output port of the interferometer (see figure 2) 
to replace the coherent vacuum state that would other- 
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FIG. 10. Graphical display of the spacetime observation vol¬ 
ume for a GW search, S = f S dt around the time of the 
OMC isolation upgrade. The period displayed here starts 
from sometime before the commissioning period until after 
all benefit-delay times have been reached. All spacetime ob¬ 
servation volumes are shifted so they are equal at the be¬ 
ginning of the commissioning period which is designated by 
the grey filled in area. Here all real spacetime observation 
volumes are plotted as solid lines while extrapolated or pre¬ 
dicted spacetime observation volumes are plotted as dashed 
lines. In black we plot quantities which correspond to the 
before-commissioning period without vetoes while red denotes 
those quantities with vetoes. The green lines plot spacetime 
observation volumes during and after commissioning (hypo¬ 
thetical and real). The benefit-delay times are displayed by 
the vertical lines where the hypothetical continued before- 
and veto-spacetime observation volume meets the after- or 
predicted-spacetime observation volume. 




FIG. 11. In the phase (P) and amplitude (A) space of the 
light field inside an interferometer we show the error ellipses 
for quantum fluctuations of a coherent vacuum state (dashed 
lines) compared to those of squeezed states (solid lines). The 
three plots show different example directions in which quan¬ 
tum fluctuations can be suppressed by squeezing. A Michel- 
son interferometer makes measurements in the phase quadra¬ 
ture of the light field so that we can represent the level of 
shot noise of the beam at the output of the interferometer 
schematically by the width of the ellipse in the P direction. In 
the region where the interferometer’s sensitivity is limited by 
shot noise the best sensitivity improvement is then achieved 
by squeezing the ellipse in the P direction as shown in (a), 
(b) shows squeezing along a direction which gives the same 
shot noise as the coherent vacuum state, i.e. where there is no 
effect from the squeezing, (c) shows squeezing along the direc¬ 
tion which increases the shot noise level in the interferometer 
output and thus decreases its sensitivity to GWs. 


wise be entering this port. If the squeezing is oriented 
correctly in relation to the electric field inside the inter¬ 
ferometer it will improve the shot noise limited sensitivity 
of the interferometer [17]; see figure 11. The orientation 
of the squeezing error ellipse with relation to the electric 
field inside the interferometer is controlled to optimise 
the improvement for the shot noise limited sensitivity as 
shown in figure 11. 

With GEO 600 leading the way in the implementation 
of squeezed light we have come across previously unseen 
noise sources. A population of transient noise events, 
“glitches”, was observed that produced broadband noise 
in the strain data channel. The noise is short duration 
and broadband above a few hundred Hz. An example of 
this can be seen in figure 13. These glitches only occur 
during periods when the squeezer is in operation. 

It was found that these glitches were due to the 
squeezer malfunctioning and producing noise that was 
being injected into the interferometer, caused by noise 
in the squeezer ellipse orientation control loop. Hence 
the error point of this loop was a viable witness chan¬ 
nel for this noise source, and was used to produce vetoes 
to remove the data contaminated by these glitches. An 
investigation correlated the low-passed time series of the 
squeezer error point to the amplitude and rate of glitches 
in the strain data channel. A threshold on the voltage of 
this channel was used to produce vetoes which removed 
the majority of the glitches, with high efficiency. These 
vetoes incur a 1.85% down-time during science opera¬ 
tions. The vetoes generated removed time intervals with¬ 
out frequency dependence as a result of the broadband 
nature of the noise contamination. 

The effect of this veto on the spectrum and cumulative 
glitch rate are shown in figure 12(a) and 12(b) respec¬ 
tively. The strain noise spectral density is not obviously 
effected—although the transient noise events are large 
in amplitude they do not happen frequently enough to 
make a significant contribution when averaging over the 
period of data used for this investigation. In the cumu¬ 
lative glitch rate histogram, the removal of these glitches 
causes the SNR threshold to be reduced by a factor of 2.5. 
In table I we see that the vetoes, by removing this noise 
source, improve the burst-like spacetime observation rate 
from 100 kpc 3 to 1200 kpc 3 . 

Although we did not obtain a full commissioning solu¬ 
tion for this problem, we can make a prediction of the per¬ 
formance of the detector with this noise source removed. 
To predict the spacetime observation rate without this 
noise source present, we use the data with our veto ap¬ 
plied but omit the 1.85% down-time it incurs. This yields 
a spacetime observation rate of 1300kpc 3 . 

Using these values of the spacetime observation rate we 
find that the predicted benefit-delay ratio from removing 
this noise source completely is 0.08. However when first 
applying the veto, the predicted benefit-delay ratio from 
removing the noise source completely is 53. This increase 
shows that the veto performs very well and in many sit¬ 
uations would make commissioning not beneficial to the 
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(a) Strain noise spectral density (b) Cumulative glitch rate 

FIG. 12. Squeezer glitches comparison plots. A period of time when glitches caused by the squeezer were happening is shown 
in black, while the same period with vetoes applied is shown in red. (a) The transient noise events are loud, but they do not 
happen frequently enough to make a significant change to the strain noise spectral density, when averaged over the period 
of data used for this investigation, (b) Both the vetoes and commissioning do well to remove the large number of high SNR 
glitches that were being produce by this noise source. The horizontal dashed black line shows the false-alarm rate threshold of 
5.6 x 1(T 4 Hz. 



Time [minutes] 


FIG. 13. Omegamap showing the transient noise events 
caused by the squeezer. 


observing capabilities of a science run. 

Without performing our cost-benefit analysis, the 
severity of this noise source as seen in figure 13 and 12(b) 
could encourage a rash decision to perform commission¬ 
ing and investigative work to attempt to locate and re¬ 
move the source of this noise. However, after applying 
our cost-benefit analysis we see that any down-time in¬ 
duced by such work would produce a deficit in our ob¬ 
servable spacetime volume which would only be compen¬ 
sated by the benefits of a fix in a period of time 53 times 
as long as the down-time. Since it was not at all clear 
how long of a down-time was needed to investigate the 
problem and commission a solution, any investigations 
which would carry on for more than a few days without 
successfully finding the source of the noise would start to 
be in danger of never being able to collect the benefits. 

In this example, our cost-benefit analysis has allowed 
us to quantitatively determine if a noise problem is oc¬ 
curring often enough to warrant commissioning. The in¬ 


frequency of the squeezer glitches and the existence of a 
good veto turned a severe problem into one that is not 
worth taking the observation time to investigate. For our 
next and last example, we describe a case where ideas 
similar to the ones presented in this article were applied 
in a decision-making process which ultimately had an im¬ 
pact to the commissioning actions. 


C. 3.5 Hz dither squeezer glitches 

A second class of glitches associated with the squeezer 
were observed during the S6E/VSR4 science run in 2011. 
Sinusoidal modulations are applied to two of the three 
BDOs to two angular degrees of freedom per optic. Each 
of these four dithers are at different frequencies. The 
purpose of the dithers is to create error signals used to 
automatically maintain the alignment of the interferom¬ 
eter beam to the output mode cleaner. A 3.5 Hz dither, 
specifically, modulates the degree of freedom with the 
largest geometric effect on the squeezer path. 

When the 3.5 Hz BDO dither was active and the 
squeezed light injection was operating, the rate of low 
amplitude glitches in the shot noise limited frequency re¬ 
gion increased, as seen in figure 16(a). Even without com¬ 
missioning work, these glitches were removed whenever 
the shutter between the squeezer and the interferometer 
was closed, and later after the 3.5 Hz dither was deacti¬ 
vated. The glitch mechanism was not entirely clear, how¬ 
ever there is a strong correlation of glitch occurrence with 
the dither period. This can be seen in figure 16(b). Most 
compellingly there was no drift in this behaviour over 
several months of operation, indicating that the glitch 
time modulation must have an origin synchronised to a 
digital clock, in this case, the real-time control system in 
which the dither sinusoid originated. These glitches do 
not show up in the Omegamaps because they were of low 
SNR, and for the same reason they also do not affect the 
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FIG. 14. Graphical display of the spacetime observation vol¬ 
ume for a GW search, S = f S dt around the time of a hy¬ 
pothetical 2 day commissioning interruption aimed at remov¬ 
ing the squeezer glitches. The period displayed here starts 
from sometime before the commissioning period until after 
all benefit-delay times have been reached. All spacetime ob¬ 
servation volumes are shifted so they are equal at the be¬ 
ginning of the commissioning period which is designated by 
the grey filled in area. Here all real spacetime observation 
volumes are plotted as solid lines while extrapolated or pre¬ 
dicted spacetime observation volumes are plotted as dashed 
lines. In black we plot quantities which correspond to the 
before-commissioning period without vetoes while red denotes 
those quantities with vetoes. The green lines plot spacetime 
observation volumes during and after commissioning (hypo¬ 
thetical). The benefit-delay times are displayed by the verti¬ 
cal lines where the hypothetical continued before- and veto¬ 
spacetime observation rate meets the predicted-spacetime ob¬ 
servation volume. 


spectrum. 

This noise source cannot be mitigated by the construc¬ 
tion of a veto as the glitches occur continuously in time 
and broadly in frequency. The only way to remove this 
glitch effect is to remove the source in the instrument 
or turn off the squeezer, an action which we will de¬ 
note as a “veto” for this investigation. Because this oc¬ 
curred during the S6E/VSR4 science run in 2011, de¬ 
tector observation time was extremely valuable and even 
without recourse to sophisticated astrophysical figures of 
merit, commissioning down-time was reserved for neces¬ 
sary maintenance and critical noise issues only. At the 
time, only the elevated glitch rate with squeezer opera¬ 
tion was known, and the dependence on the 3.5 Hz dither 
did not become clear until the last weeks of the science 
data run. Early in the science run, the question arose 
as to whether to disable the squeezer, to incisively in¬ 
vestigate the source of the glitches, or to allow squeezer 
operations to continue and accept the additional low SNR 
transients. 

Fortunately during the science run, the ability to carry 
out studies similar to those described in this paper was 
available. Elements of the coherent wave burst all-sky 


burst search analysis [18] were processing the coincident 
GEO600-Virgo network data on a weekly basis. These 
analyses provided a figure-of-merit for the non-stationary 
elements of the network noise similar to p in the in¬ 
tegrated fixed false-alarm probability range defined by 
equation 1 and equation 2. A directed analysis of the 
3.5 Hz dither squeezer glitches which compared time in¬ 
tervals with the squeezer shutter open and closed showed 
that the effect of these glitches on the network fixed false- 
alarm probability SNR threshold was less than a few per¬ 
cent. In the single interferometer study shown in figure 
15, we obtain similar results. Since the squeezing at that 
time was acting to decrease the noise spectral density 
of GEO 600 by approximately 20%, we decided that no 
intervention involving the sacrifice of science data opera¬ 
tions was warranted. We see the effect of this in table I; 
there was a ~ 25% loss in the burst-like spacetime obser¬ 
vation rate when the squeezer was off which clearly shows 
that deactivating the squeezer is counter productive. 

To predict the performance of the detector without 
this noise source we combined the spectra of the detec¬ 
tor with squeezing, and the fixed false-alarm probability 
SNR threshold found without squeezing. This gives us an 
idea of how the detector would perform with squeezing 
on and without the noise source. This gave us a predicted 
burst-like spacetime observation rate of 1700 kpc 3 . 

Using these values we find a benefit-delay ratio when 
comparing a time using squeezing to our predicted per¬ 
formance without the noise source of 33; this means we 
would have had to find the problem and fix it within 
two days for us to be able to collect the benefits within 
the two months left in the run. Thus performing in¬ 
vestigations or attempting to commissioning this issue 
would risk inducing a deficit in potentially observed GW 
source which could not be recovered before the end of 
the S6E/VSR4 science run. The decision to not commis¬ 
sion was based on results from the coherent wave burst 
analysis performed at the time of the noise, but utilising 
the full benefit-delay analysis we obtain a quantitative 
measure of the risk that was involved in commissioning. 


VI. SUMMARY 

As we have seen in the course of characterisation in¬ 
vestigations of GEO 600, different detector noise sources 
require different mitigation strategies to most effectively 
remove the noise induced reduction of the spacetime ob¬ 
servation rate. When these noise sources occur during a 
science data run, we compare the benefit-delay time when 
we first apply the veto solution to the remaining time in 
the observing period. If the commissioning down-time 
plus the benefit-delay time is smaller than the remaining 
time in the observing period, then the commissioning 
should be performed to maximise the spacetime obser¬ 
vation of the science data run. Otherwise it is better to 
apply any available vetoes and then perform the commis¬ 
sioning after the observing period has ended. Using this 
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FIG. 15. 3.5 Hz squeezer glitches comparison plots, (a) A period before hypothetical commissioning with the noise source 
present is shown in black, and a period with out veto, removing the squeezer, is shown in red. (b) We see that this noise source 
has no notable effect on the rate of triggers in the SNR range frequency band, from 500 Hz-4 kHz. The horizontal dashed black 
line shows the false-alarm rate threshold of 5.6 X 10 -4 Hz. 



(a) Cumulative glitch rate for triggers above 1 kHz 


(b) Histogram of Omega trigger time modulo 2 seconds 


FIG. 16. Triggers in a frequency band from 1 kHz-4 kHz, where this noise source can be observed, (a) We see that in the 
frequency band from 1 kHz-4 kHz the effect of these low SNR glitches is more clearly visible. The horizontal dashed black line 
shows the false-alarm rate threshold of 5.6 x 10 -4 Hz. (b) Histogram of trigger time modulo 2 seconds. This noise source can 
clearly be seen in the before period, with the squeezer active, but is not observed in the before with veto period when the 
squeezer was turned off. 



Investigation 

Range type 
[frequency 
band] 

Burst-lil 

Before 

(x 10 5 kpc 3 ) 

ce spacetim 

Vetoed 

Before 

(xlO 5 kpc 3 ) 

Le observati 

Predicted 

(x 10 5 kpc 3 ) 

on rates 

After 

(x 10 5 kpc 3 ) 

Beam Direction 
Optics (BDO) 

GRB 

[100 Hz-500 Hz] 

6 

45 

300 

100 

Output mode 
cleaner (OMC) 

GRB 

[100 Hz-500 Hz] 

470 

750 

850 

1240 

Squeezer glitches 

SN 

[500 Hz-4 kHz] 

0.001 

0.012 

0.013 

- 

3.5 Hz dither 

squeezer glitches 

SN 

[500 Hz-4 kHz] 

0.016 

0.012+ 

0.017 

- 


TABLE I. Burst-like spacetime observation rates for each of the investigations given in section V. These spacetime observation 
rates are calculated from the astrophysically motivated fixed false-alarm probability ranges as described in section II and are 
displayed for a period before commissioning, the same period before commissioning with vetoes applied, the predicted results 
of a commissioning operation, and for a period after a commissioning operation was completed, 
f The veto time for the 3.5 Hz dither squeezer glitches was to turn off the squeezer. 
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Investigation 

Range type 
[frequency 
band] 

Commissioning 

down-time 

(days) 

no vetoes, 
predicted 

Benefit- 

no vetoes, 

actual 

delay ratios 
with vetoes, 
predicted 

with vetoes, 

actual 

Beam Direction 
Optics (BDO) 

GRB 

[100 Hz-500 Hz] 

20.8 

0.02 

0.06 

0.2 

0.75 

Output mode 
cleaner (OMC) 

GRB 

[100 Hz-500 Hz] 

1.9 

1.3 

0.6 

8 

1.5 

Squeezer glitches 

SN 

[500 Hz-4 kHz] 

- 

0.08 

- 

53 

- 

3.5 Hz dither 

squeezer glitches 

SN 

[500 Hz-4 kHz] 

- 

33 

- 

2.6 f 

- 


TABLE II. Commissioning down-times and benefit-delay ratios for each investigation in section V. We calculate the benefit- 
delay ratios, as described in section II equation 5, for scenarios where commissioning down-time is used to bring about the 
improvements to the spacetime observation rates shown in table I. Here we look at multiple cases for each investigation where 
we either take the spacetime observation rate calculated without vetoes or with vetoes for the era before commissioning, and 
either the predicted or actual spacetime observation rate for the era after commissioning. The total benefit-delay time can be 
calculated by multiplying the ratios with the actual commissioning down-times, 
f The veto time for the 3.5 Hz dither squeezer glitches was to turn off the squeezer. 



FIG. 17. Graphical display of the spacetime observation vol¬ 
ume for a GW search, S = f S dt around the time of a hypo¬ 
thetical 2 day commissioning interruption aimed at removing 
the 3.5 Hz dither squeezer glitches. The period displayed here 
starts from sometime before the commissioning period until 
after all benefit-delay times have been reached. All spacetime 
observation volumes are shifted so they are equal at the be¬ 
ginning of the commissioning period which is designated by 
the grey filled in area. Here all real spacetime observation 
volumes are plotted as solid lines while extrapolated or pre¬ 
dicted spacetime observation volumes are plotted as dashed 
lines. In black we plot quantities which correspond to the 
before-commissioning period without vetoes while red denotes 
those quantities with vetoes. The green lines plot spacetime 
observation volumes during and after commissioning (hypo¬ 
thetical). The benefit-delay times are displayed by the verti¬ 
cal lines where the hypothetical continued before- and veto¬ 
spacetime observation rate meets the predicted-spacetime ob¬ 
servation volume. 


cost-benefit figure-of-merit we can quantitatively com¬ 
pare any available vetoes with the predicted improvement 
from a commissioning solution, to best inform decisions 
on detector commissioning. 

In section V A we investigated noise in the alignment 
between the interferometer and the OMC. A veto which 
removed the frequency bands contaminated by this noise, 
and a commissioning solution to remove the source of 
this noise were available. In retrospect, for the BDO 
suspension and OMC isolation upgrades we predicted 
benefit-delay ratios which predicted short benefit-delay 
times compared to the science data run duration, and 
can verify that this was achieved in the end, though the 
actual ratios varied from the predicted by factors of a 
few. 

In section VB we discuss a situation in which the 
glitches are short in duration and broadband, and an 
auxiliary channel is available from which a time segment 
veto can effectively be produced. Initially the noise ori¬ 
gin was not well known, and there was no commissioning 
solution available to remove this noise, but we could pre¬ 
dict the performance of the detector without this noise 
source present. The benefit-delay ratio found compar¬ 
ing the veto to our predicted performance without this 
noise source is 53, indicating that the benefit-delay time 
would have been 53 times the down-time for any com¬ 
missioning. Any investigative work for a commissioning 
solution that caused down-time would also require a large 
benefit-delay time, so it is favourable to perform non in¬ 
vasive investigations or work at the same time as other 
commissioning activities. Therefore, applying the veto 
for this noise source was the best solution available. 

In section V C we see that the 3.5 Hz dither squeezer 
glitches were not significantly affecting the detector’s as- 
trophysical sensitivity. Work similar to this investiga- 
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tion, but instead using a network figure-of-merit which 
included the Virgo detector, was performed when this 
noise source was observed and was key in making the 
decision to maintain operations with the squeezer in op¬ 
eration, although it was the known source of this noise. 
The cost-benefit analysis shows that the improvement in 
the spectrum, due to the squeezer, has a larger effect than 
the increase in the false-alarm rate threshold, due to the 
3.5 Hz dither squeezer glitches, on the burst-like space- 
time observation rate. Therefore the best solution was 
continue to operate with the squeezer as it was, rather 
than disabling the squeezer or interrupting the observa¬ 
tions to investigate further. 

In this article, we have successfully developed, refined, 
and applied a cost-benefit strategy to approaching com¬ 
missioning decisions aimed at solving noise sources in 
GEO 600 observed during and around the time of the 
S6E/VSR4 data science run. Since the development of 
this strategy happened only during and after the science 
run, the application in its final state, presented here, was 
only done after the fact. For the upcoming era of ad¬ 
vanced GW detectors, however, it is important to apply 
this strategy to commissioning decisions, as they arise, 
before and during scientific observing runs. Accurate es¬ 
timation of the commissioning time is necessary for the 
successful application of this strategy, which in our ex¬ 
perience has been available from the experience of the 
commissioners themselves. 

There is one further step which needs to be taken be¬ 
fore application in the advanced detector network can 
become a reality. All of the spacetime observation rates 
reported here have been calculated based on a fixed false- 
alarm probability range for a single detector, externally 
triggered, unmodelled GW burst search. In the advanced 
detector era, however, the flagship GW search which is 
expected to make a detection is a multiple detector net¬ 
work search for modelled waveforms from compact bi¬ 
nary coalescences (CBCs). It is important that a fixed 
false-alarm probability range is developed for this search. 
This is a non-trivial extension of the ranges used in this 
article, however, much progress has already been made 
toward reducing the computational delay of generating 
false-alarm probabilities for CBC triggers from a network 
of two or three interferometers [ ]. With this hurdle 

overcome, it should be simple to construct a fixed false- 
alarm probability range for CBC searches on a network 
of interferometers modelled after the ranges presented in 

[i ]• 

Once this is done, the cost-benefit strategy to com¬ 
missioning decisions developed here can be carried out 
for commissioning decisions on individual interferometers 
within a network. The effects of these decisions will be 
tied into the entire network by the network range figure- 
of-merit. Of course with a network there will be more op¬ 
tions for commissioning because the question may arise 
whether or not it makes sense to take down more than 
one interferometer at a time. However, the strategy in 
principle remains the same. 


In the advanced detector era, characterisation groups 
can and must evaluate the sensitivity cost of noise phe¬ 
nomena as they appear, develop vetoes, and compare the 
relative effectiveness of vetoes and commissioning solu¬ 
tions on the network of interferometers. This exciting 
role provides a vital contribution in the preparation and 
operation of the first advanced detector observation runs. 
With the detectors in pursuit of design sensitivity and 
beset by technical noises, this strategy will help to max¬ 
imise the total spacetime observation volume and bring 
forward the first direct detection of GWs. 
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